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Abstract 

We apply Pade summation to the /3(A) function for the quartic Higgs coupling 
A in the standard electroweak model. We use the (5 function calculated to five 
loops in the minimal subtraction scheme to demonstrate the improvement re- 
sulting from the summation, and then apply the method to the more physical 
on-mass-shell renormalization scheme where (3 is known to three loops. We 
conclude that the OMS (3 function and the running coupling A(^) are reliably 
known over the range of energies and Higgs-boson masses of current interest. 
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I. INTRODUCTION 



In the standard model of elementary particle physics, the SU(2) x U(l) symmetry is 
spontaneously broken to a residual U(1) EM , generating mass for the W ± and Z gauge 
bosons and the matter fields. A possible cause for the symmetry breaking is the presence 
of an additional scalar field, the Higgs field. Although there is as yet no experimental 
evidence for the expected Higgs boson, we can still explore the implications of this symmetry 
breaking mechanism using radiative corrections to standard-model processes. For example, 
the condition that perturbative calculations be reliable provides a theoretical upper bound 
on the mass Mh of a weakly interacting Higgs boson [|T|-|TT||. 

In a recent paper ||TT[] , Nierste and Riesselmann analyzed one- and two-scale processes 
involving the Higgs field with a particular emphasis on the running of the quartic Higgs 
coupling A(/i). They assessed the reliability of perturbation theory using two criteria: the 
relative difference of physical quantities calculated in different renormalization schemes; and 
the dependence of A on the renormalization scale [i. If perturbation theory is to be reliable, 
the choices of the renormalization scheme and scale should not be important for physical 
quantities. To determine A(/z) in their analysis, Nierste and Riesselmann integrated the 
renormalization group equation using the three-loop f3 function, and solved the resulting 
equation for A iteratively using four different approximation schemes. The solutions differed 
significantly for large values of the coupling or mass scale, and determined one constraint on 
M H in a perturbative theory. This uncertainty in A(/i) carries over to physical quantities such 
as scattering amplitudes and again affects the ranges of Mh and fi over which perturbative 
calculations are reliable. 

We show here that is is possible to explore the regime of large coupling without the 
ambiguities that arise from the direct iterative solution for the coupling. We approach the 
problem by emphasizing the (3 function /3(A), and show that it can be determined reliably 
to rather large values of A by using Pade approximates |12|,|13J to sum the perturbation 
series for /3(A) |14|]. Integration of the renormalization group equation then gives an implicit 
equation for A(/z) that can be inverted numerically. The results can be used to study the the 
validity of perturbation theory for scattering amplitudes in the region of large Higgs-boson 
masses and high energies where the running coupling A(/i) is the natural renormalization 
group expansion parameter. We will not pursue those applications here as a number of 
authors [^- [TT| , p!5[] have considered them in detail. 

We first investigate the Pade approach in Sec. |H| using the the results for /3(A) in the 
minimal subtraction (MS) renormalization scheme for which the perturbation series for (3 
is known to five loops [pl|[l7]]. After establishing the effectiveness of the Pade approach, 
we apply it in Sec. |nj to the more physical on-mass-shell (OMS) renormalization scheme 
where (3 is only known to three loops |pT| , [18| . We find that Pade summation of the series 
apparently gives a reliable result for (3{\) for quite large values of the coupling, A < 10, 
and conclude, after inversion of the renormalization group expression, that A(/z) is known 
reliably in the OMS scheme for < 4 TeV for Mh < 800 GeV. The region in which A(/z) 
is known well extends to very large mass scales if Mh is sufficiently small, for example, to 
10 17 GeV for M H < 155 GeV. 
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II. PADE SUMMATION OF THE (3 FUNCTION 



A. Preliminary considerations 

In the following, we deal with the quartic Higgs-boson coupling A defined at tree level 
in terms of Mh and the electroweak vacuum expectation value v = 246 GeV or the Fermi 
coupling Gp by A = M^ 2 /2v 2 = GfMh 2 / y2. We will work in the interesting limit of large 
Higgs-boson masses, corresponding to the limit of large quartic couplings, and neglect the 
effects of couplings with fermions. 

The running coupling A(/i) is defined as the solution of the renormalization group equa- 
tion 

/-^ = m m 

at the energy scale fi. The function /3(A) is given in perturbation theory as a power series 
in A, 



\2 / A \ n 



The coefficients f3 n are renormalization-scheme dependent beyond two loops. They are 



known through three loops in the on-mass-shell renormalization scheme |TT, 

OMS : fa = 24, fa = -312, (3 2 = 4238.23, (4) 
and to five loops in the minimal subtraction scheme |16|J17|| , 



MS : (3 = 24, /3i = -312, (3 2 = 12022.7 

p 3 = -690759, /3 4 = 4.91261 x 10 7 . (5) 

Alternatively, the coefficients B n are given by 

OMS : Bq — 1, B x = -0.082323, B 2 = 0.0070816, (6) 

MS : B = I, B l = -0.082323, B 2 = 0.020089 

£ 3 = -0.0073090, Bi = 0.0032917. (7) 

To determine the running coupling, one must integrate the renormalization group equa- 
tion, Eq. (HD, and solve the implicit equation 

In— = / -wr-y 8 
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This equation determines A(/i) in terms of the initial and final mass scales /io and /i and 
the initial value of the coupling at the scale /z , defined as A = A(/i ). Different, typically 
iterative, methods of solution lead to different results for A(/i), with the differences increasing 
for large vales of M H or A and for ji^> /i . Since the j3 function is only known to finite 
order, the only constraint on this standard approach is that the different solutions satisfy the 
renormalization group equation, Eq. (§), to that order. However, the resulting ambiguities 
for large values of Mh can compromise tests of the reliability of perturbation theory, and 
the determination of limits on Mh in a weakly interacting theory. It is therefore useful to 
approach the problem differently, and concentrate on the /3 function itself. If /3(A) is known 
accurately for some range of A, the integral in Eq. ([8]) will also be accurately determined, 
and the equation can be inverted numerically to find A(/i) in that region. 



B. Pade summation and /3(A) 



Pade approximates fli~2| , |i~3f give a very useful way of summing or extrapolating series for 
which only a finite number of terms are known. The [N, M] Pade approximate for a function 
f(z) defined by a truncated power series 

m 

f(z) = Y / c^+0(z m+1 ) (9) 

3=0 



is a ratio of two polynomials, 



E 



N 



P[N, M](z) = ^T° : 6 ° = 1 ' N + M = m. (10) 

Z^ n =0 

The coefficients determined uniquely by the requirement that the series expansion 

of P[N, M](z) agree term-by-term with the series for f(z) through terms of order z m . 

The sequence of Pade approximates P[N, M] is known to converge to f(z) as N, M — > oo 
with N — M fixed for large classes of functions |l2l , |l"3| , but the approximates can also give 



useful and rapidly convergent asymptotic approximations for finite N and M even if the 



sequence and the original series for f(z) do not converge [13|. 

In the present case, the function in question is /3(A), known perturbatively to orders A 4 
and A 6 , that is, to three and five loops, in the OMS and MS renormalization schemes, re- 
spectively. The perturbation series for /3 is not expected to converge, but a Pade summation 
of the series may still be useful for A not too large. Because the perturbative expansion 
of /3(A) starts at order A 2 , we will extract the leading power explicitly, redefine the Pade 
coefficients, and define the [N, M] approximate for the n-loop (3 function as [] 

m , r n A 2 l + aiA + a-2A 2 H V a N \ N „ T _ . . 

mM\=^ t + _ „ , N + M = n-l. (11) 



x Pade summation of was considered by Yang and Ni [14|, but without applications to the 
present problem. Those authors did not extract the overall factor A 2 , so use a different labeling of 
the approximates, and miss the diagonal approximates used here. 
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Note that the approximates (3 [n — 1, 0] are just the perturbation series for (3 carried to n 
loops. 

The series for /3(A) defined by Eq. (|3]) are alternating series in which the ratios of co- 
efficients B n+ i/B n change only slowly in either OMS or MS renormalization in the range 
in which the S's are known. This suggests that the diagonal approximates (3 [N, N] with 
M = N or the subdiagonal approximates with M = N + 1 may be particularly effective 
in estimating the series. In the case of OMS renormalization, the /3's are known only to 
three loops, so M + N < 2. The possible choices are then (3 [1, 1] or (3 [0, 2] if we use all 
the three-loop information, or (3 [0, 1] if the perturbation series is truncated at two loops. 
f3 [2, 0] and (3 [1, 0] are just the three- and two-loop perturbation series. In the case of MS 
renormalization, (3 is known to five loops, M + N < 4, and we will consider the approximates 
(3 [1, 2] at the four-loop level, and (3 [2, 2] at five loops, keeping M = N or M = N+1. 
The additional five-loop approximates (3 [1, 3], (3 [3, 1], and (3 [0, 4] are members of sequences 
two or more steps off the diagonal. These are not expected to converge as rapidly as the 
sequences we consider. The coefficients aj, bj for these approximates are given in appendix 
A. 



C. Tests of Pade summation using MS renormalization 

The fact that the perturbation series for (3 is known to five loops gives us the opportunity 
to test the Pade summation procedure using known results. Having established its reliability, 



we will the apply the method in Sec. Ill to the more physical OMS renormalization scheme 



in which the connection between A and M# is known. 



1. Convergence of the Pade sequence 

Based upon the general convergence properties of Pade approximates and the alternating 
character of the series at hand, we expect the sequence (3 [1, 1], (3 [1, 2] and /3 [2, 2] to 
converge as we progress from three to five loops. We plot these approximates in Fig. ^ to 
demonstrate that convergence. The convergence of the Pade sequence is, in fact, relatively 
fast. For low values of A there is excellent agreement. Even for A = 10, (3 [1, 1] and f3[l, 2] 
differ by < 10 % with the diagonal five-loop approximate (3 [2, 2] lying roughly halfway 
between the other two. We interpret the agreement and the pattern of convergence as 
strong evidence for the effectiveness of the (3 [N, N] sequence in summing the series for 
(3(X), and conclude that it is unlikely that (3 would be found to differ significantly from 
(3 [2, 2] in the region shown if higher- loop contributions were calculated. 

In Fig. ^ we look at the problem from the point of view of the purely perturbative 
approach, and show the sequence of the N-loop perturbation series (3 [N — 1, 0] for f3. This 
is not a sequence in which A" and M increase together with the difference N — M fixed, so the 
standard results on Pade convergence do not apply. The convergence of the sequence is very 
slow as shown in the figure, with large differences between successive terms already present 
for A ~ 3. For comparison, we also show the three- and five- loop diagonal approximates 
f3[l, 1] and (3 [2, 2]. These forms interpolate the perturbative sequence very well, eliminating 
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the dominance of the last term in the series for A large. Since /3[1, 1], fi [2, 2], and the four- 
loop approximate fi [1, 2] differ from each other by less than 5% for A < 10, all are effective 
in extrapolating the perturbation series. We conclude, in particular, that the three-loop 
approximate /3[1, 1] already gives a reliable extrapolation for /3(A), with uncertainties of 
only a few percent, out to A ~ 10, far beyond the range in which the five-loop perturbation 
series is reliable. 



2. Estimates of unknown coefficients 

Pade approximates often converge to the limit function faster than the power series used 
to construct them. In that case, the terms in the expansion of a Pade approximate beyond 
the matched order may give reasonable estimates for the unknown higher-order coefficients 
in the power series. As a simple test of this expectation in the present case, we can expand 
the three- and four- loop approximates (3 [I, 1] and (3[1, 2] to one order higher in A than the 
finite power series used to construct them, and compare the new coefficient with the known 
four- and five-loop results. Thus, the expansion 

fi [1, l] = -^(3 \l + B x \ + B 2 X 2 + (B 2 2 /B 1 )X 3 + (B 3 2 /B 2 )\ 4 + ■ ■ ■} (12) 
gives the estimates 

B' 3 = B 2 /B 1 , B'± = B\jB\, (13) 

for the four- and five-loop coefficients B 3 and B4, results equivalent to 

/3' 3 = fill fix = -463286, fi' A = fiUfil = 1.785 x 10 7 . (14) 

The actual four- and five-loop results are 

fi 3 = -690759, fi A = 4.913 x 10 7 . (15) 

The estimates of fi 3 and fi& from the three-loop coefficients are therefore about 0.67 and 0.36 
of the actual coefficients. 

In the case of fi [1, 2], we can estimate only /3 4 , with the result 

B' A = - (Bl - 2B X B 2 B Z + B 3 2 ) j (B\ - B 2 ) . (16) 



This estimate gives fi' 4 = 3.48 x 10 7 , and a ratio fi'J fi± = 0.71. 

The estimates for the first missing terms in the perturbation series are too small in both 
of the cases considered. We can understand this result qualitatively as resulting from the 
averaging of an alternating series by the approximates, with the corresponding tendency to 
avoid large higher coefficients in the expansion. We will use this observation below. 

The effects of incorrect estimates of B% on the approximate fi[l, 2] are shown in Fig. 
|3|. In these calculations, we have taken B 3 as five- and ten times the estimated value, and 
calculated fi [1, 2] using the new value as input. The result is a < 10% change in fi for 
A < 10 despite the very large values of the new coefficient. 
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D. The running coupling A(/x) in MS renormalization 



The effect of the uncertainty in /3(A) on the running of A(/i) can be studied by integrating 
the renormalization group equation, Eq. (p]), and solving numerically for A as a function of 
its initial value Ao and the ratio of energy scales ////^o-Q We have done this calculation using 
the approximates j3 [1, 1] and (3 [1, 2], choosing initial values A = 1, 3, 5. The results are 
shown in Fig. |j. The result for the optimum five-loop approximate, (3 [2, 2], lies near the 
center of the shaded regions in that figure, as would be expected from the comparison of 
the approximates in Fig. |l|. We believe the estimated range of uncertainty is quite generous 
given the rapid convergence of the sequence shown there toward (3 [2, 2]. 

The range of uncertainty in A(/i) at fixed /V/^o is quite small for Ao = 1, 3 over the entire 
range shown, /x//i < 6. The uncertainty is larger for A = 5, roughly 16%, at /x///o — 3, but 
even then the boundary curves differ from the curve for (3 [2, 2] by < 8%. 

The rather small effect of uncertainties in (3 on A(/i) can be understood rather simply. 
The renormalization group equation involves 1/(3 rather than (3. The prefactor A 2 in the 
Pade expression in Eq. ([□]) leads to a rapid decrease in the integrand, and the value of the 
integral is determined mainly by the region near Ao, the lower endpoint of the integration. 
For Ao small, (3 is well determined in the most important region, and the uncertainty in the 
integral is small. The uncertainty in the integral, hence the uncertainty in A(/i), becomes 
large only for renormalization group evolution away from a large starting value for Ao- 



III. APPLICATIONS: RANGES OF RELIABILITY OF /5(A) AND \(p) IN OMS 

RENORMALIZATION 

A. Pade approximates for f3(\)oMS 

Having tested the use of Pade approximates in the MS scheme, we consider the impli- 
cations of Pade summation for the OMS scheme. The most significant difference is the 
limited order, three loops, to which the perturbation series for (3 is known. We are therefore 
restricted to two approximates that use the full information available, the diagonal approx- 
imate (3 [1, 1] and the subdiagonal approximate (3[0, 2]. We can also use (3[Q, 1] at the 
two-loop level. Based upon the convergence of the Pade sequence demonstrated for the MS 
scheme, and the apparent reduction in the size of the coefficients in the OMS scheme^, we 
will assume that these approximates again provide an accurate estimate for the (3 function, 
with the diagonal approximate probably the most reliable. 



2 In the case of MS renormalization, A is connected only indirectly to the physical pole mass of 
the Higgs boson, so we cannot state the results in terms of \i and Mh without using a separate 
calculation of the self-energy function. 

3 The known value of /?2 in the OMS scheme is smaller than that in the MS scheme by roughly a 
factor of three [ [Tl]]T8[ ]. Nierste and Riesselmann [11] have found similar reductions in the coefficients 



in the expansion of physical amplitudes. We assume that the reductions in the size of the coefficients 
persist at higher orders. 
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To determine the range of A for which the (3 function is reliable, we first considered the 
differences among the three-loop functions (3[1, 1] and (3 [0, 2], and the two-loop function 
(3 [0, 1]. These approximates can barely be distinguished over the range of A shown in Fig. |5| 
with the scale used there, so only (3\1, 1] is shown. This agreement is the result of the 
nearly geometric growth of the first coefficients in the perturbation series. The three-loop 
approximates j3 [1, 1] and (3 [0, 2] continue to agree well to much larger values of A. While 
one is tempted on this basis to conclude that the OMS (3 function is reliably known for 
A < 10, the range of current interest, the geometric character of the low-order perturbation 
series may well be accidental. We have therefore attempted to estimate a wider range of 
uncertainty in the f3 function in a different way by supposing, in agreement with the results 
of the MS analysis, that the coefficient B' 3 estimated by expanding (3 [1, 1] is too small, and 
constructing a new "four-loop" approximate (3[1, 2] using a greatly increased value of B 3 . 
The result obtained using B 3 = 5B' 3 is shown in Fig. |5|. The change in the extrapolation 
of the perturbation series is quite small, with a difference of less than 2% between f3[l, 1] 
and (3 [1, 1] for A < 10. We also show the perturbation series for the (3 function, (3 [2, 0], in 
Fig. [5] for comparison. 



B. The running coupling A(//) 

In the OMS renormalization scheme, the parameter A is defined by the relation A = 
GfMh 2 /\/2 to all orders in perturbation theory ||^0|. We will choose the starting value 
A of the running coupling A(/i) to have this value. What remains to be decided is the 
energy scale fiQ at which this relation should be taken to hold. The natural energy scale 
would appear to be /i = M H . However, other choices have been made. Thus, in an early 
investigation, Sirlin and Zucchini |IJ calculated the one-loop corrections to the four-point 
Higgs-boson scattering amplitude and defined the parameters in the theory so that large 
electromagnetic effects appear only in such standard relations as that between Gf and the 
muon decay rate. With this definition, the high mass limit of the four-point function gives 

m 



h{n) = A 



1 + 



Ao 
16vr 2 



24 In 



M; 



^— + 25 - 3\/37r 



H 



(17) 



The logarithm in the expression above is just that which appears in the expansion of the 
one-loop expression for A(/z), 



A(/i) = A c 




(18) 



for (3q = 24 and fi = Mh- The ambiguity in the choice of [i is in the treatment of the 
remaining constants in Eq. (|T7|) . These have been incorporated in the running coupling 
by some authors [p^ ,p|,p| by redefining fi as /^o = Mffexp[(— 25 + 3\/37r)/24]. However, 
the constants do not appear naturally in the expression for the four-point function at two 
loops [[J. It is probably most reasonable, therefore, to treat them as separate "radiative 
corrections" and write h(u) to one loop as /i(/i) = A(/z)[l + <5], with A(/i) the one-loop running 
coupling defined above, and 5 incorporating the remaining scale-independent corrections. 



S 



This question has been studied in more detail by Nierste and Riesselmann |11]], who 
showed that the convergence of the perturbation series was improved for several physical 
amplitudes by adopting the natural scale /i = M H instead of the choice noted above. They 
note, furthermore, that in order to cancel large logarithmic terms in the perturbative result 
when one considers two-scale physical processes such as scattering, the scale /i must be 



related to the energy scale of the interaction by /x = |TTJ. We will follow Nierste and 
Riesselmann and make the definite, physically motivated choices /Jq = Mh and fi = y/s in 
the following analysis. This specification amounts as already noted to a definite specification 
of the "radiative corrections" in perturbatively calculated amplitudes once the couplings are 
expressed in terms of A(/x). 

With Ao and fiQ specified, and the range of reliability of the (3 function established, 
it is straightforward to integrate the renormalization group equation and invert the result 
numerically to obtain A(/f). The uncertainty in A(/x) can be specified in terms of that in (3. 



With this procedure, no further uncertainties such as those illustrated in [II] are introduced. 
It is not necessary, for example, to obtain the solution of the renormalization group equation 
as a series in Ao, in an iterative approximation. We note in this connection that the "naive" 



and "consistent" forms for A(/x) given in |TT| correspond respectively to the approximates 



(3[N, 0], the perturbation series for (3, and (3 [0, N], the series obtained by expanding (3. 
Neither sequence is expected to converge well with increasing N. 

Our results for A(/i) are shown in Fig. |6] for M H = 500 and 800 GeV and /i < /■* < 4 
TeV. We find for Mh = 500 GeV that all Pade approximates, including the perturbation 
series, agree very for /j, < 5 TeV, a region in which Ao < 5. The residual uncertainty in A(/x) 
is small enough not to affect perturbative results for physical processes. 

Different Pade approximates also give very similar extrapolations for A(/j) for Mh = 800 
GeV, even when the predicted value of (3% is changed by a large factor. The only significant 
deviation involves the perturbation series (3 [2, 0] which we do not believe is reliable on the 
basis of our earlier investigation. Even if we restrict the range of A in which we take (3 as 
reliable to A < 10 as in Fig. [5], the result for A(/i) remains reliable for /i < 2 TeV, a value well 
into the energy region of interest for experiments at the Large Hadron Collider at CERN. 

The rapid growth of A(/i) for the perturbative approximate (3 [2, 0] in Fig. ^] is the result 
of an Landau pole at fi — 2339 GeV. A pole can appear in A(/i) if the integral of 1/(3 
converges for A — > oo, with the position of the pole in fi determined by the condition 

ln^=hm/ A ^. (19) 

Hq a^oo J\ [3 

No pole can actually appear when the integration is restricted to the finite range of A in 
which f3 is known reliably, but the likely presence of a pole would be indicated by very rapid 
growth of A(/x) with increasing /i in that region. In the present case, there is no reason 
to expect the perturbation series (3 [2, 0] to be accurate for A large. The results in Fig. |5| 
indicate, in fact, that the perturbative approximation begins to fail badly for A ~ 5, while 
the starting point for the evolution of A(/i) shown in Fig. ^| is at A = 5.3 for Mh = 800 GeV. 
The remaining approximates do not lead to poles in the region shown. 
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C. Conclusions 



We have shown that Pade summation of the (3 function improves the reliability of (3 and 
the running quartic Higgs coupling. The method gives a best estimate for (3, and removes 
much of the uncertainty associated with different determinations of A(/i) at the three- loop 
level [|TTJ. 

We have tested the Pade method using the (3 function in the MS renormalization scheme, 
where f3 is known to five loops in the perturbation expansion. The test results suggest rapid 
convergence of the diagonal and subdiagonal Pade sequences. Our applications are to the 
more physical OMS renormalization scheme, where the first scheme- dependent coefficient 
in the OMS expansion is significantly smaller than in the MS expansion. This more rapid 
apparent convergence is reflected in the excellent agreement among the leading Pade ap- 
proximates for Poms even for rather large values of the first unknown coefficient, (3%. 

Application of the Pade method to the three-loop results for the OMS (3 function leads to 
a running coupling that appears to be quite reliable in the region studied, namely a Higgs- 
boson mass M H < 800 GeV and a mass scale ji < 4 TeV. We conclude that uncertainties in 
A(/i) are not an important source of uncertainty in perturbative results for physical scattering 
and decay amplitudes in this interesting region. 
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APPENDIX A: PADE APPROXIMATES 

The coefficients of the Pade forms used in our analysis are given below. We will state 
the results in terms of the coefficients B n = (/? n //3 )/(167r 2 ) n . 

At two loops, N + M — 1 and we have only the truncated perturbation series f3 [1, 0] 
and the approximate (3 [0, 1] with 

/? [0, 1] : b l = -B l (Al) 

At three loops, N + M = 2 and we have the new approximants (3 [1, 1] and (3 [0, 2}. The 
coefficients are given by: 

P[l,l]: a 1 = (B 2 1 -B 2 )/B 1 , (A2) 
h = -B 2 /B x , 

£[0,2]: b 1 = -B 1 , (A3) 
b 2 = B\ — B 2 . 

At five loops, N + M = 4 and we will consider the new approximants f3 [3, 1], (3 [2, 2], 
(3 [1, 3] and (3 [0, 4]. The coefficients are given by: 
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/? [3, 1] : a 1 = (B 1 B 3 -B 4 )/B 3 , (A4) 
a 2 = (B 2 B 3 — BiB 4 )/B 3 , 
a 3 = (Bl - B 2 B 4 )/B 3 , 
61 = -B 4 /B 3 , 



(3 [2, 2] : ai = (5i£ 2 2 - B\B 3 + 5x5 4 - B 2 B 3 )/A 22 , (A5) 
02 = (-B2 ~~ ^B\B 2 B 3 + B\Bk + .83 — B 2 B 4 ) / A 22 , 
6i = (5i# 4 - B 2 B 3 )/A 22 , 
b 2 = (I? 2 — -B2-B 4 )/A 2 2, 
A 2 2 = B\ — BiB 3 , 



(3 [1,3]: ai = (B? - 35^2 + 2S1S3 + ^ - 5 4 )Mi 3 , (A6) 
h = {—B\B 2 + B\ + - S 4 )/A 13 , 
b 2 = [B\B\ — B 2 B 3 — B\B 3 + B 1 B 4 )/A 13 , 
h = (2BiB 2 B 3 — B 2 — B\ + B 2 B± — B x B^)j A i3 , 
^-13 = B\ — 1B X B 2 + S 3 . 

We will also consider the approximates 2], the subdiagonal approximate for the 
four-loop expansion. The coefficients in this case are: 

(3 [1, 2] : ai = (5? - 25^2 + 5 3 )/( J B 1 2 - fl 2 ), (A7) 
b 1 = (B 3 -B 1 B 2 )/(B 2 1 -B 2 ), 
b 2 = (Bl-B 1 B 3 )/{B 2 1 -B 2 ). 



APPENDIX B: ANALYTIC RESULTS 



The Pade approximates we have used are all integrable analytically. We will give only 
the results needed in our investigation of the OMS renormalization scheme: 



5 2 - 25 2 -Bi + 2S 2 A 

H — = arctan =, (B2) 

/45 2 - 5 2 JAB 2 - B\ 



( R2 _ d \ 

+£? 1 ln(l + 1 2 Aj, (B3) 
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These expressions are to be equated to ((3 /l6ir 2 ) ln(/i//! ). 



12 



REFERENCES 



[1] D.A. Dicus and V.S. Mathur, Phys. Rev. D 7, 3111 (1973). 

[2] B.W. Lee, C. Quigg, and H.B. Thacker, Phys. Rev. Lett. 38, 883 (1977); Phys. Rev. D 
16, 1519 (1977). 

[3] M.S. Chanowitz and M.K. Gaillard, Nucl. Phys. B261, 379 (1985). 
[4] W. Marciano, G. Valencia, and S. Willenbrock, Phys. Rev. D 40, 1725 (1989). 
[5] S. Dawson and S. Willenbrock, Phys. Rev. Lett. 62, 1232 (1989); Phys. Rev. D 40, 2880 
(1989). 

[6] G. Passarino, Phys. Lett. 156B, 231 (1985); Nucl. Phys. B343, 31 (1990). 
[7] C.E. Vayonakis, Europhys. Lett. 12, 29 (1990). 

[8] L. Durand, J.M. Johnson, and J.L. Lopez, Phys. Rev. Lett. 64, 1215 (1990); Phys. Rev. 
D 45, 3112 (1992). 

[9] L. Durand, P. Maher and K. Riesselmann, Phys. Rev. D 48, 1084 (1993). 
[10] L. Durand, B.A. Kniehl and K. Riesselmann, Phys. Rev. Lett. 72, 2534 (1994); (E) 74, 

1699 (1995); Phys. Rev. D 51, 5007 (1995). 
[11] U. Nierste and K. Riesselmann, Phys. Rev. D 53, 6638 (1996). 
[12] G. Baker, Jr., Essentials of Pade Approximants, Academic Press, 1975. 
[13] C. Bender and S. Orszag, Advanced Mathematical Methods for Scientists and Engineers, 

McGraw Hill, 1978. 
[14] J.-F. Yang and G.-J. Ni, Comm. Theor. Phys. 22, 207 (1994). 
[15] A. Ghinculov, Phys. Lett. B 337, 137 (1994); (E) 346, 426 (1994). 
[16] A. A. Vladimirov, D.I. Kazakov and O.V. Tarasov, Sov. Phys. JETP 50, 521 (1979). 
[17] H. Kleinert, J. Neu, V. Schulte-Frohlinde, K.G. Chetyrkin and S.A. Larin, Preprint 

hep-th/950320 (1995). 
[18] M. Liischer and P. Weisz, Phys. Lett. B 212, 472 (1988). 
[19] A. Sirlin and R. Zucchini, Nuclear Physics B 266, 389 (1986). 
[20] L. Durand and K. Riesselmann, Phys. Rev. D 55, 1533-1547 (1997). 



13 



FIGURES 



FIG. 1. A comparison of the sequence of five-loop diagonal and subdiagonal Pade approximates 
for /3(A) in the MS renormalization scheme. Note that alternate approximates are too large or too 
small, and that the sequence converges rapidly with the final result presumably in the band between 
(3[1, 2] and f3[2, 2]. 

FIG. 2. Demonstration of the slow convergence of successive perturbative approximations to 
the (3 function toward the diagonal Pade approximate (3 [2, 2] for MS renormalization. 

FIG. 3. Plots of the Pade approximate (3 [1, 2] to the MS [3 function using the actual value of 
the four-loop coefficient (3^ and values five and ten times the estimate obtained from the three-loop 
approximate f3[l, 1]. 

FIG. 4. Plots showing the running of A(/x) as a function of the ratio of scales n/ no for different 
initial choices of Ao in the MS renormalization scheme. The differences between the curves obtained 
using the Pade approximates f3 [1, 1] and (3 [1, 2] corresponding to three- and four-loop summations 
of (3 indicates the range of uncertainty in the result. The curves for the five-loop (3 function lie 
near the center of the band of uncertainty. 

FIG. 5. Plots of the two- and three-loop Pade approximates (3 [0, 1] and (3 [1, 1] for (3 in the 
OMS scheme. The function f3 [1, 2] obtained using a coefficient f3s five times as large as that 
estimated from (3 [1, 1] is shown to indicate a range of uncertainty. The three- loop perturbation 
series (3 [2, 0] is shown for comparison. 

FIG. 6. Plots showing the running of A(/Lt) in the OMS renormalization scheme as a function 
of the mass scale n for different initial choices of the Higgs-boson mass Mh ■ The differences 
between the curves obtained using the three-loop Pade approximate (3[1, 1] and the function (3[1, 2] 
obtained using a coefficient (3% five times as large as that estimated from (3 [1, 1] is shown to indicate 
a range of uncertainty. The perturbative result for (3 given by (3 [2, 0] has a Landau pole at n = 2.3 
TeV for Mh = 800 GeV, but is not reliable and is included only to illustrate the effects of a nearby 
pole. 
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